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We study the exciton condensate (EC) in a bilayer two-dimension-electron-gas (2DEG) adjacent 
to a type-II superconductor thin film with an array of pinned vortex lattices. By applying con- 
tinuum low energy theory and carrying numerical simulations of lattice model within mean-field 
approximation, we find that if the order parameter of EC has a vortex profile, there are exact zero 
modes and associated rational fractional charge for zero pseudospin potential (/i) and average chem- 
ical potential (h): /i=0 and /i=0; while for and /i=0, intervalley mixing splits the zero energy 
levels, and the system exhibits irrational fractional axial charge. 
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Introduction Charge fractionalization and fractional 
statistics [1] in two-dimensional systems have attracted 
much attention due to their potential applications in 
topological quantum computation [2-4]. One representa- 
tive example is the well known i^=| fractional quantum 
Hall (FQH) system [Hj. The excitations of the iy=^ FQH 
system carry fractional charge and exhibit non-Abelian 
fractional statistics. The pivotal features of FQH system 
are strong correlations and broken time-reversal symme- 
try (TRS) as a result of the electrons' strong coulomb 
interaction and strong external magnetic field. However, 
some recently proposed models [6, 7] have revealed that 
the strong correlations and broken TRS are not neces- 
sary for fractionalization. Moreover, these models indi- 
cate that the nontrivial topological configurations play 
a key role in the fractionalization. The common fea- 
tures of these models are that the TRS is not broken, 
and the low energy excitations can be described by the 
Dirac field coupling with a topologically nontrivial scalar 
field (vortex profiles, for instance) [8|. Therefore, the 
energy spectrum's symmetry and index theorem of the 
Dirac operator protect the zero modes and the associ- 
ated fractionalization [9]. In another interesting model 
[lO] different from aforementioned ones, the fractional- 
ization of the excitations is due to the response to the 
half-magnetic-flux-quanta defect of the system in the in- 
teger quantum hall (IQH) regime. This model is more 
likely to be fabricated compared to aforementioned sys- 
tems. 

We start with a bilayer two-dimension-electron-gas 
(2DEG), as illustrated schematically in Fig. 1(a). Each 
layer is adjacent to a film of type-H superconductor sup- 
plying the quantization of flux in unit of {^o=h/e) 
by the Abrikosov vortex square lattices. The spin- 
polarized electrons in the bilayer are almost localized 
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FIG. 1: (Color online) (a) Schematic structure of exciton 
condensate system. The external gates induce the symmet- 
ric electron and hole carriers, (b) Our bilayer square lattice 
model. 



near the Abrikosov vortex cores. At half filling (namely, 
the number of electrons is one half of the total lattice 
sites), the bilayer system has a magnetic filling factor of 
two, and the magnetic flux through one square lattice 
cell is TT. In addition, the two layers are separated by a 
spacer which is usually a dielectric barrier (e.g., Si02). 
The chemical potential in each layer can be adjusted inde- 
pendently by the respective gate voltage. The carrier in 
each layer has electron or hole nature determined by the 
respective band structure. In particular, the inverse ex- 
ternal gate voltages in two layers bias the charge balance 
of two layers and result in the perfect particle-hole sym- 
metry. So the electrons in one layer and holes in another 
layer bound together and form exciton condensate (EC) 
due to the interlayer coulomb interaction. Our system 
is favorable for EC because of the strong magnetic field 
and low electron density as suggested by experiments of 
the coupled quantum wehs [ll| [12| . 

In this paper, we construct a tight-binding model anal- 
ogous to bilayer graphene [13] and solve the Bogoliubov- 
de Gennes (BdG) equation self-consistently in real space 
to determine the amplitude of EC order parameter (OP) 
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and study the system's topological properties originating 
from OP with vortex configuration. The low energy ex- 
citations of the system present novel features due to the 
vortex configuration in OP. We show that when the bi- 
ased voltage is zero, two degenerate fermion zero modes 
(each one for one valley) emerge. Furthermore, for the 
vortex with vorticity (n=l) at half filling, each layer has 
fractional charge or — (depending on whether the 
two zero modes are occupied or not), though the total 
and difference charge of bilayer are the integer numbers. 
When the biased voltage is not zero, the two fermion 
zero modes slightly and symmetrically split around zero 
energy due to the intervalley mixing. At half filling, the 
axial charge of the whole system is fractional just as what 
happens in bilayer graphene studied in Ref. (13). 

Lattice model and exciton condensation In the 
following, we present a model to describe the EC in our 
system at near zero temperature. The essential terms 
include intralayer hopping between the nearest neigh- 
bor sites. We ignore the intralayer carriers' Coulomb 
interactions because their effects only renormalize the en- 
ergy bands. We consider inter layer Coulomb repulsion as 
an effective short-range repulsion only between interlayer 
electrons at the same planar sites. Furthermore, the ex- 
ternal biased voltage is included. Since the electron spin 
is polarized along the direction of the magnetic field, we 
omit spin index. There is no direct hopping between the 
two layers. 

The Hamiltonian of our system can be expressed as 
follows 



ia <i,j>a 



(1) 



Here Qq, annihilates an electron at site i of layer-a 
(a=l,2), ni, 



tot ^icx^ioi^ 



and tijcx is intralayer nearest- 
neighbor hopping with tiji=tij2 for the electron-electron 
bilayer and tiji=—tij2 for electron- hole bilayer. We con- 
sider the electron-hole bilayer and set tiji=—tij2=t=l. 
jioi is the chemical potential of layer-a. The mag- 
netic field is included through the Peierls phase fac- 
tors, Oija = ^ r^"^ Aq,-^!. We choose to work in Lan- 
dau gauge Aq, = ^<I>o(— ^, 0), where we have set the lattice 
space unity. The mean-field order parameter of the EC 
is Ai=U (c^Qi) , and the mean-field Hamiltonian is 



ia <ij>a 

-^(A,4c,2 + /l.c.) + ^^|A,|^ (2) 

i i 

With the help of the following Bogoliubov transformation 
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the above Hamiltonian (j2j) can be diagonalized by solving 
the following BdG equation 



E 



-A*Sij hij^2 



AiSij 



En 



(4) 



where /iij,i(2) = — Mi(2)<^ii- The self-consistent equa- 
tion of the OP is 



Ai = tanh( 



En . 

2kBT' 



(5) 



We solve the set of BdG equations self-consistently 
via exact diagonalization in real space. The system 
size of 24x24 is used in the calculation and the con- 
vergence criterion of is set to 10~^ in unit of the 
intralayer nearest-neighbor hopping. We consider the 
near-zero temperature and set ksT to be 10~^. For 
convenience of discussion, we define the average chem- 
ical potential and pseudospin polarization potential as 
/i=(/ii+/i2)/2 and /i=(/ii— /i2)/2. We find that if h=0 
(we set /i2=— Mi=— m), the exciton order parameter is 
uniform: Ai=A=m, where m is real. Our calcula- 
tions show that as h decreases from zero, the phases of 
system undergo BCS-like phase, ID modulated Larkin- 
Ovchinnikov-like (LO) phase [1^, 2D modulated LO-like 
phase and finally OP zero. 

In the following, we restrict ourselves to the case of 
/i2=— Mi=— In the Landau gauge that we choose, the 
system can be viewed as simplified lattices as depicted 
in Fig. 1(b). The magnetic unit cell in layer a in- 
cludes two sites (/, n)a and (/, n-\-l)a denoted by Aa and 
Ba. The wave functions are of the form of spinor field 

V^(r) = [cBi(rz,n+i),CAi(r^,n),CB2(r^,n+i),CA2(r^,n)]^. In 
the momentum space with the reduced Brillouin zone 
BZ={ky—7T/2<kx<ky-\-7T/2, \ky\ <7r}, the Hamiltonian 

can be written as i^MF=X]k '^^(^)'^^^(^)'^(^)+^o 
with Eo= Nqw? /U and 

HmfO^) = -j^os (g) /+2 cos{kx)as 53 - 2 cos{ky)as si 



(Jl /, 



(6) 



where (si) (i=l,2,3) is Pauli matrix acting in sublat- 
tice (pseudospin) space and / is unitary matrix. With- 
out loss of generality, we preserve the constant phase x 
of OP. This Hamiltonian possesses the particle-hole sym- 
metry 



nuMF{\<-)n = -uiip{-'i<-) 



(7) 



with Q=—(j2^I and Q^=l. As a result of this symmetry, 
it is easy to find 



(8) 



It means that if i/j is an eigenvector of Hmf with eigen- 
value E, then ftip'^ is guaranteed to be an eigenvector of 
Hmf with eigenvalue —E. So, all nonzero eigenvalues of 
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Hm f come in ±E pairs and the energy spectrum is sym- 
metric about zero energy. For uniform order parameter 
m, the energy spectrum is of the form 



^(k) = ±2y (^cos2 A:^ + cos2 ky ± iif + m?. (9) 

The energy spectrum symmetry is clear in this case. Note 
that no in-gap eigenvalues exist as m^O in uniform case. 
Zero modes and charge fractionalization We now 

investigate the low-energy behaviors of our system. In 
the low energy limit, the properties of the system are 
dominated by the excitations around the two inequiv- 
alent valleys K± = (^^,^^) where o^o^ kx\j^^j^^=^Px 
and cos ky\^^^^=±.py. Then, the linearized mean-field 
Hamiltonian for these excitations are 
with l-LMF+=-'iI^cF2^MF-'iI®cF2=l-L, 

U = -/iCr3(g)/+2pa,a3(8)53-2p2^Cr3(8)Si-me~'^'^3^^Cri(g)/. 

(10) 

Here, p=— zV is the momentum operator. Since the val- 
leys are decoupled in the low-energy theory, we first study 
a single valley. 

Around the valleys, the Hamiltonian 1-i possesses a 
symmetric property 

S^(K±+p)5 = ^*(K±+p), (11) 

where S=i<j2(8)S2. This symmetry is unbroken only 
around each valley. The low-energy excitations obey 
this symmetry around each valley. The case of /i=0 is 
of particular interest. In this situation, the symmetry 
of our system gets enhanced. It is easy to check that 
Ql-L9=-n, 6^=1, with 9=cr30S2. This symmetry has 
important consequences on the number of zero modes. 

Now suppose A is also taken to vary with position 
and consider OP with a symmetric vortex of vorticity n, 
which means A(r, 0)=m{r)e^^^ with m(r) a real function 
of r and the constant phase x omitted by the special 
gauge choice. In order to get explicitly analytic solutions, 
we use the well defined vortex profile 

In quantum limit, m(r)=mo, where mo is a positive real 
constant. The zero- mode solutions satisfy HV^o=0- Eq. 
(pT]) means HS?/^*=-SH*V^*=-£^S?/;*. For zero modes, 
we impose S?/^q=?/;o. If S?/;q=— ?/^o, then the transforma- 
tion i/jQ^i^l^o brings it back to EiljQ=ipo. So, in the zero- 
energy subspace, we can suppose ipo = [f,g^ig^^—if*]^. 
Hence, we can get two independent equations for zero- 
mode solutions in quantum limit 

(-/i - 2id,)f + 2idyg - zmoe^"^^* = 

(-/i + 2id,)g + 2idyf + imoe^^V* = 0. (13) 

In order to get the symmetric solutions in the presence 
of the symmetric vortex profile, we introduce new wave 



functions ( and ^ defined by the following equations. 
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Then, satisfy the fohowing equations in polar coor- 
dinate, 

e'^dr + ir-^deX - mo/2e'"''C* - «/^/2? = 
e-'^'idr - ir-'^dg)i + mo/2e"^r - ifJ'/K = 0, (15) 

which are almost the same as that in Ref. (13). The 
analogous solutions can be gotten as 

C = e-™<"'/2 Jp(Mr/2)e*f'' and ^ = -ie-'""''/^ J,(/xr/2)e'''^ 

(16) 

which are exponentially decayed Bessal functions 
and p=q—l=^^^^. So, the single valued solu- 
tions survive only when n is odd. Hence, one 
zero-mode solution ipQ^\r) can be obtained as 

4'HO=[4i]i(0,^^'ii(r),^^i]2(^).^o'i2(0]^=[/,5,i5*,-i/1^ 
together with Eqs. (p!4|) and (p!6|) . From the relation 

l-LMF-\-=—iI'S)(J2l-LMF-iI^o-2^ another zero- mode so- 
lution is i^f\r)=[ig,-if, f* ,g*Y . Note that the 
particular symmetry (defined by the transformation O) 
for /i=0 ensures the \n\ zero modes, unlike the single 
zero mode for odd vortex configuration for ji 7^0 [18]. 

With the help of Eqs. (pj) . (fT5)) and (fTG]) . it 
is interesting to see that the low energy behav- 
ior of our model is similar to that of bilayer 
graphene. The analogy is not surprising, because 
after the transformation /iMF(k)=5'^HMF(k)S' with 
5'=exp(i7r/(8)5i / 4) exp(i7r/(g)53 / 4) exp(i7r/(g)5i /2) , our 
mean-field Hamiltonian has the form 

/iMF(k) = -7o(7i2cos/c^+722cosA:2^+me~*^^^+7o75/i). 

(17) 

Here in the Weyl representation, the Dirac ma- 
trices have the forms: 'y^=io'2^s^^ 7o=cri(8)/, 
7^=— i7o7i7273=cr3(8)/. The low-energy form of 
this Hamiltonian at the two independent valleys is 
nearly identical to that for bilayer graphene. 

A localized zero mode in a gapped system with 
particle- hole symmetry usually bounds a fractional 
charge ±e/2 [lej. Note that the Hamiltonian 

of our model in low-energy limit is of the Dirac type 
with first-order differential operators. Therefore, when 
/i=0, and the OP has the form of a symmetric vortex 
of vorticity n, i.e., Eq. (fT2]) . there are \n\ independent 
zero modes existing for each Dirac valley indicated by 
the index theorem [17]. When jn^O^ the special symme- 
try related to transformation © disappears and there is 
a single zero mode (given by by Eqs. ([16]) and (p!4|) ) for 
odd vorticity as argued in [18]. Other zero modes for 
|n| >1 associated with one Dirac valley must split from 
zero energy and form in-gap bound states which partially 
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(a) (b) 




FIG. 2: (Color online) (a) Charge density distribution (with 
the uniform background charge density distribution sub- 
tracted) for n=l, /i=0 and at half filling case with two zero 
modes occupied. The charge (the integral of the charge den- 
sity) is equal to |e. Here we give one layer case, and the other 
layer has the same distribution, (b) The axial charge density 
distribution with n=l and 0.3. In all cases, we set the 
lattice system size as A/'=32x32x2, h—0, and U—A. 



preserve zero modes properties just as pointed out in Ref 
(9). Of course in general the inter- valley mixing would 
split the pairs of zero modes. In this case, only axial 
charge fractionalization may occur, which is in analogy 
with that in Ref. (13). The case of /i=0 is an exception 
due to the symmetry. The two zero modes belonging to 
different valleys are degenerate and do not split in spite 
of the existence of inter-valley mixing. Hence, each layer 
of the vortex core bounds a fractional charge ±e/2 when 
the two zero modes are occupied or not. Our numerical 
calculation shown in the next section verifies the above 
picture. 

Numerics We perform exact diagonalization of the 
mean- field Hamiltonian of our system with 32x32 sites 
per layer for symmetric vortex profile of the OP in 
quantum limit with the amplitude of OP calculated 
from the self-consistent solution of the BdG equations. 
In Fig. 2 we show the charge density distribution 
[Fig. 2(a)] and the axial charge [Fig. 2(b)] defined by 
^Qz^=Xli[(^^i)zy ~ (^^i)o]5 with subscript u{0) referring 
to the system with 1(0) vortex and Sni=ni2—nii the dif- 
ference in the occupation number of the two layers. 

The numeric results show the rational charge faction- 
alization for /i=0 and irrational fractional axial charge 
for /iT^O. Here the rational fractional charge is exactly 
±e/2, independent of U. The dependence of SQ^ on /i 
and U in Fig. 3 (a) shows that the axial charge may be 
continuously tuned by changing the external parameters. 

The splitting of the zero energy levels due to inter- 
valley interaction is shown in Fig. 3(b). The ab- 
sence of splitting of zero energy levels for /i=0 can 
also be understood from a perturbative analysis. It 
can be seen that the energy splitting has the form 
SEoc\moJ e-'^^+'^-^'^Rele-'^ {f^g^)]dr\. One can see 
that ^^=0 with the help of solutions of Eqs. (14) and 
(16). Therefore, the inter- valley coupling does not split 
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FIG. 3: (Color online) (a) The axial charge bound to a vor- 
tex profile as a function of U and — /x. (b) Energy splitting 
[E{N/2+l) - E{N/2)] changes as a function of U and -/i. In 
both cases, /i=0 and the vortex profile has vorticity n=l. 



the degenerate zero modes. The existence of the ex- 
act zero modes ensures the charge fractionalization (per 
layer) as seen in Fig. 2 (a). 

Discussions We propose a system of bilayer 2DEG 
adjacent to a type-II superconductor thin film for EC. 
The vortex profile for OP can be generated by defects or 
additional nontrivial texture in magnetic field p^]. Fur- 
thermore, in our bilayer system, if a vacancy lies on a 
site in one layer and meanwhile an interstitial induced by 
impurities occurs at the underlying site in another layer, 
the nontrivial exciton OP with vortex profile is generated 
due to its minimal coupling to the nonzero axial gauge 
field emerged from the vacancy and interstitial relative 
to the magnetic fiux lattices background. Moveover, the 
axial fractional charge bounding a magnetic flux quanta 
can form charge-flux composite particles which behave as 
Abelian anyons with the exchange phase j=7r6Qu under 
the standard argument [20] . Note that the above studies 
have ignored the intralayer next-nearest-neighbor hop- 
ping_ti which may be important in the similar systems 
[2l|, [22I , where the fractional charge is related to topolog- 
ically nontrivial gapless edge states. However, the charge 
factionalization discussed in the present paper is different 
from that case, since the inter-layer interaction destroys 
the gapless edge states. Besides, unlike the crude crys- 
talloid, the influence of ti here can be negligibly small 
since the lattice distance of the Abrikosov vortex square 
array can be adjusted very large. 

The bilayer EC is analogous to the superconductor by 
the particle-hole transformation. That is why the BdG 
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equations were used. Furthermore, many phenomena in 
superconductor have counterparts in bilayer EC, e.g., the 
vacancy induced suppression of the amphtude of the OP 

In summary, we have studied EC in a bilayer system 
which is hkely to be fabricated in laboratory. In the 
presence of topological non-trivial OP profiles (vortex 
profiles), we have shown that the system exhibits ratio- 



nal/irrational fractional (axial) charge. 
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